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Abstract. We present a combined radio, optical, and X-ray study of the nearby LINER galaxy NGC 1052. Data from a short 
(2.3ksec) CHANDRA observation of NGC 1052 reveal the presence of various jet-related X-ray emitting regions, a bright 
compact core and unresolved knots in the jet structure as well as an extended emitting region inside the galaxy well aligned 
with the radio synchrotron jet-emission. The spectrum of the extended X-ray emission can best be fitted with a thermal model 
with kT = (0.4 - 0.5) keV, while the compact core exhibits a very flat spectrum, best approximated by an absorbed power-law 
with Nu = (0.6 - 0.8) x 10"^ cm"-. We compare the radio structure to an optical "structure map" from a Hubble Space Telescope 
(HST) observation and find a good positional correlation between the radio jet and the optical emission cone. Bright, compact 
knots in the jet structure are visible in all three frequency bands whose spectrum is inconsistent with synchrotron emission. 
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NGC 1052 is a nearby' elliptical galaxy which harbors a low- 
luminosity active galactic nucleus (LLAGN) in its very center 
(Li-iooGHz = 4.4 X 10'*"ergs-i; Wrobel H^S^ . It hosts a two- 
sided radio jet emanating from the nucleus and reaching out to 
kiloparsec-scales which is, however, still fully enclosed within 
the stellar body of the optical galaxy. In the optical, the spec- 
trum of NGC 1052 is characterized by strong forbidden lines 
from low-ionization states which has made NGC 1052 the pro- 
totypical LINER (low-ionization nuclear emission line region; 
Heckman 1980) galaxy. As for LINERs in general, it has long 
been argued whether these low-ionization Unes in NGC 1052 
are excited by a central photo-ionizing source (e.g., Gabel 
et al. EUOO ) or if shock heating is the dominant mechanism 
(e.g., Sugai <& Malkan EOOOI I. While there is overwhelming ev- 
idence for the presence of an active galactic nucleus (AGN) in 
NGC 1052, the role of shocks in this galaxy is still unclear The 
improved angular resolution in X-rays offered by CHANDRA 
makes it possible for the first time to image the distribution of 
X-ray emission on the same scales as accomplished by con- 
nected radio interferometers, e.g., MERLIN (Multi-Element 
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Radio-Linked Interferometer Network). Disentangling the con- 
tributions of compact (i.e., <1 arcsec) nuclear and extended 
(i.e., >1 arcsec) X-ray emitting regions to the total amount of 
X-ray emission of NGC 1052 can serve as an important tool 
to study the interaction between the radio jet plasma and the 
ambient interstellar medium. 

NGC 1052 has been observed by all major X-ray mis- 
sions of the pie-CHANDRA era, like Einstein (Mc Dowell 
[r994j . ASCA and ROSAT (Weaver et al. .1999j . and Beppo 
Sax (Guainazzi & Antonei mi999> . For these X-ray missions, 
NGC 1052 appeared as a point-like X-ray source. The X-ray 
spectrum of NGC 1052 is extremely flat, a finding that led to 
the proposal of an advection dominated accretion flow (ADAF) 
as the origin of the observed X-ray emission (Guainazzi et 
al. 2000 1. To model the AGN X-ray spectrum above E ^ 
2keV absorbing column densities in excess of lO^-' cm"^ have 
been discussed, supporting the idea of a high density obscur- 
ing torus, comparable to column densities found in other AGN 
(e.g., Mahzia et al l 19971 Risahti et al. l2002> . Independent ev- 
idence for the existence of an obscuring torus at the center of 
NGC 1052 is obtained from Very Long Baseline Interferometry 
(VLBI) observations in the radio regime: on parsec-scales 
NGC 1052 exhibits a twin jet structure with a prominent emis- 
sion gap between both jets (see e.g., Kadler et al. 2003b I. 
The inner part of the western jet shows a strongly inverted 
radio spectrum, which was first discovered by Kellermann et 
al. (.1999 J (see also Kameno et al. [200 Ij . The cm-wavelength 
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0.2-1.0 keV 



Fig. 1. Raw CHANDRA image of NGC 1052 showing the spatial dis- 
tribution of photon events in the energy range (0.2 - 1.0) keV. 



spectral index in this central region is larger than 2.5, exceeding 
the theoretical limit for synchrotron self-absorption. 

Combined studies of the core region of NGC 1052 in the 
radio and X-ray regime are of essential importance for con- 
straining the physical properties of the parsec-scale radio jet 
and the obscuring torus as well as to determine the nature of 
the nuclear X-ray emission. In this paper we present a com- 
bined radio, optical, and X-ray study of the jet-related emis- 
sion in NGC 1052 on arcsecond scales. In particular, we focus 
on the soft X-ray excess in the source-spectrum below E - 
2keV. This soft component was identified first by Weaver et al. 
(fT999^) based on ROSAT PSPC data. CHANDRAs superior an- 
gular resolution makes it possible to present evidence that this 
soft excess emission is associated with the well known radio 
jet. 

In Sect.|2lwe present the CHANDRA, MERLIN, and HST 
data as well as their reduction. In Sect. |3l we discuss the 
arcsecond-scale morphology of NGC 1052 in the radio, opti- 
cal, and X-ray regime and the correlations between the differ- 
ent wave bands. In Sect.0]we derive models for the nuclear and 
extended X-ray emission and Sect. |5] summarizes our conclu- 
sions. 

2. Observations and data reduction 

2.1. CHANDRA data 

CHANDRA observed NGC 1052 on August 29/30, 2000l 
During the 2342 s observation, the Advanced CCD Imaging 
Spectrometer (ACIS) Chip S3 was in the focus of the High 



^ The CHANDRA data were taken from the public archive 
jhttp : // 'cxc ■ harvard . edu/cda/chaser . html ( and analyzed us- 
ing standard methods within the software package ciao 3.0 using the 
calibration data base caldb 2.2.5 version. 



Resolution MiiTor Assembly (HRMA). The ACIS-S3 detector 
offers high angular resolution as well as information on the X- 
ray source spectrum because of its intrinsic energy resolution 
and has a higher sensitivity in the soft X-ray energy regime 
below £■ < 1 keV then the front-side illuminated detectors of 
CHANDRA. 

The nucleus of NGC 1052 is the brightest X-ray source 
within the field of interest. The count rate of 0. 12 cts s ' is suf- 
ficiently high to affect the measured AGN X-ray spectrum by 
the pile-up effect. The pile-up effect changes the shape of the 
measured X-ray spectrum as well as the measured count rate, 
because during a single read-out period of the chip multiple 
X-ray photons may be detected within a single pixel. Due to 
insufficient time resolution their combined signal is registered 
as a single photon event. The presented CHANDRA observa- 
tion was performed using the standard timed exposure mode. 
Depending on the detailed shape of the AGN X-ray spectrum 
up to 27% of the available data can be affected by the pile-up 
effect. In principle, this introduces a bias which results in spec- 
tral hardening and mimics lower count rates. 

The spatial distribution of photon events on the ACIS-S3 
chip during the observation of NGC 1052 is shown in Fig.^for 
the energy range (0.2 - l)keV. Even from the small number 
of counts during this snapshot observation the distribution of 
soft X-rays clearly differs from a point spread function. Using 
the csMOOTH program which is part of the ciao software we also 
produced an adaptively smoothed map of the field of interest 
(see Fig. 13 for which we reduced the angular resolution of the 
data from 0.5 arcsec to 4 arcsec and set a minimum significance 
threshold to 3cr. 

2.2. MERLIN data 

A MERLIN observation of NGC 1052 at 1.4 GHz was per- 
formed on November 22, 1995. The data from this experi- 
ment have been obtained from the public archive^ and ana- 
lyzed applying standard methods using the program difmap. 
Two different maps of the brightness distribution of NGC 1052 
at 1 .4 GHz were produced. First, a strong (u, v)-taper was used 
to map the extended emission resulting in a restoring beam of 
(1.5 X 1.1) arcsec at a position angle (RA.) of -31°. Second, a 
pure naturally weighted image was produced yielding a restor- 
ing beam of (0.5 x 0.3) arcsec at a P.A. of 27°. 

2.3. HST data 

For this paper, we re-examined the archival HST data of 
NGC 1052 previously published by Pogge et al. (2000). The 
archival data are three narrow-band F658N filter images ac- 
quired with the WFPC2 in the PCI detector (pixels scale ~ 
0.046 arcsec) with integration times of 200, 800, and 900 s, 
respectively. Cosmic rays were cleaned from the individual 
images using a version of the L.A. Cosmic Laplacian edge- 
detection algorithm described by van Dokkum (200l|l. The 
cleaned images were added together after registration using a 
field star in the image to remove the inter- image offsets. To en- 
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Fig. 2. Three color CHANDRA image 
of the jet-associated X-ray emission in 
NGC 1052 with the 1.4 GHz MERLIN 
radio image overlaid in contours. The 
CHANDRA image has been smoothed to 
a resolution of 4" to increase sensitivity. 
Red represents photons between 0. 1 keV 
and 1 ke V, green between 1 keV and 2 keV 
and blue (2 - 3)keV. The clipping level 
is 3cr. The MERLIN image was restored 
with a circular beam of 1.7 arcsec FWHM. 
Contours are shown for (0.5, 1.5, 2.5, 
3.5, 4.5, 5.5, 6.5) x 1 mjy/beam. The jet- 
associated emission is dominated by soft X- 
rays below 1 keV, whereas the AGN appears 
white, representing a rather flat X-ray spec- 
trum. 



hance the faint emission features noted by Pogge et al. ('2000V 
we created a "structure map" of the combined, cleaned images 
following the technique described by Pogge & Martini (p002 1. 
We used a model point spread function for the F658N filter 
and PCI camera generated using the TinyTim package (Krist 
& Hook 1999 1. Structure mapping is an improved method of 
image contrast enhancement (compared to traditional unsharp 
masking). In brief, the technique suppresses the large-scale 
starlight distribution, enhancing faint emission and absorption 
(e.g., dust extinction) features in the image. 

3. The multi-waveband jet structure 

3.1. Jet-associated X-ray emission in NGC 1052 

The CHANDRA image of NGC 1052 in the energy range 
(0.3 - 3.0) keV is shown in Fig.|2l Strong X-ray emission from 
the AGN is seen in this image, as well as diffuse extended emis- 
sion well aligned with the radio jet, whose MERLIN image is 
superimposed in contours. In Sect. 14.1 1 and Sect. 14.21 we dis- 
cuss the X-ray spectra in detail and present evidence that the 
soft excess component found previously in the ROSAT spectra 
of the nucleus (Weaver et al. .1999j is associated with this soft, 
extended jet-associated emission. 

To the first order, the radio and the X-ray jet are aligned 
on arcsecond scales and their extent appears to be the same. 



However, the intensity of synchrotron emission in the radio 
lobes appears to be anti-correlated with the soft X-ray intensity 
distribution in so far as the radio hotspots on both sides cor- 
respond to regions on the ACIS-S3 chip which did not detect 
any photon events. East and west of the nucleus two relatively 
bright X-ray emitting regions coincide with emission knots in 
the radio regime. Bright optical emission knots are also present 
in these areas (Pogge et al. i2000t . 

3.2. Radio emission on arcsecond scales 

The large-scale radio structure of NGC 1052 visible in the ta- 
pered MERLIN image reveals some differences to the 1980 
VLA image of Wrobel iT984b (see Fig.|3}. The VLA, which 
was operated in its A configuration during the 1980 observa- 
tions, was more sensitive to extended structures, which par- 
tially have been resolved out by MERLIN. The core has var- 
ied, with a flux density of ~0.74Jy in December 1980 and 
~1.01Jy in November 1995. The western hot spot has in- 
creased in flux density. A knot in the western jet (labeled as 
B in the MERLIN image) was not visible as a local maxi- 
mum in the 1980 VLA image. By aligning the core positions at 
both epochs, we found that the eastern knot (labeled as A), has 
moved about 0.5 arcsec inwards. However, an examination of 
the pure naturally weighted MERLIN image (see Fig.E} reveals 
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that knot A is composed of two isolated sub-components (A 1 
and A 2) separated by the same distance of ~0.5arcsec. This 
suggests that the inner sub-component might have increased 
in flux density between 1980 and 1995, causing an apparent 
shift of the blended knot structure in Fig.|2 This offers a more 
plausible explanation than an inward motion with a velocity of 
~lOc. 
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Fig. 3. Comparison between the VLA image of NGC 1052 observed 
in 1980 (Wrobel [T^M and the MERLIN image from 1995, both at 
1.4 GHz. Both maps have been convolved with the same (circular) 
restoring beam of 1.7arcsec FWHM. Contours at (-1, 1, 2, 3, 4, 5, 
6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 800) x 0.5 mjy are given in both 
maps (the VLA map also shows negative contours at (-4, -3, -2) 
X 0.5 mJy). The hot spots (labeled as HI and H2) do not show any 
positional changes. The emission peak of the eastern knot (A) is about 
0.5 arcsec closer to the core in 1995. 



3.3. Optical structure 

The structure map of NGC 1052 is shown in Fig.0]with the ra- 
dio contours of the untapered, full-resolution MERLIN image 
superimposed. The jet-counterjet structure is aligned with the 
optical emission cone. The dark band perpendicular to the ra- 
dio jet might be an artifact of the image processing. Two optical 
emission knots are located at the edges of the two radio sub- 
components of knot A. The optical emission knot in the west 
coincides roughly with a weak (~ lo") radio feature while the 
stronger radio knot 1 .5 arcsec further out has no corresponding 
bright optical counterpart. The origin of the optical emission 
remains unclear since there is no continuum image to subtract 
from the Ha filter image of NGC 1052. The conical morphol- 
ogy in this LINER 1 .9 galaxy is remarkably similar to the struc- 
ture of the narrow-emission-line region typically observed in 
Seyfert2 galaxies (e.g., Falcke et al. .l998j . This suggests that 



the optical emission cone is due to line emission rather than 
continuum emission. The optical flux density of the two eastern 
knots of (68 + l)juJy exceeds the power-law extrapolation from 
the radio to the X-ray regime by almost three orders of magni- 
tude. This is a strong argument against the synchrotron emis- 
sion process, but compatible with relatively strong line emis- 
sion exceeding the contimuum emission of the knots in this 
narrow band. 

4. X-ray imaging spectroscopy 

Using CHANDRA'S high angular resolution, it is possible to 
obtain separate spectra of the nucleus and the diffuse jet- 
associated X-ray emission. We selected an annulus centered at 
the nucleus position which excludes the nucleus itself but in- 
cludes the whole area of the diffuse emission. To avoid any 
contamination by the nuclear X-ray emission, a diameter of 
3 arcsec for the inner annulus ring was used. The nuclear X- 
ray spectrum was extracted up to a diameter of 2 arcsec. The 
extracted X-ray data of both regions were corrected for the 
unrelated X-ray background emission, using a region located 
close by without significant point sources. We also applied a 
correction for the quantum efficience degeneracy of the ACIS 
detector. 

Fig.lsjreveals the much softer spectrum of the diffuse X-ray 
emission compared to the nuclear spectrum. The bulk emission 
originates below E < 2 keV, while the X-ray spectrum of the 
nucleus has an additional hard X-ray component. Remarkable 
is the soft X-ray emission below E < 2 keV towards the nu- 
cleus. ASCA, as well as the ROSAT PSPC, could not separate 
the nucleus and the jet spatially, but a soft excess X-ray emis- 
sion had already been detected in the X-ray spectrum from the 
PSPC data. 

4.1. The X-ray spectrum of the extended emission 

To constrain the emission process of the diffuse extended emis- 
sion, we fitted different models using xspec (version 1 1.2): the 
RAYMOND model corresponding to a Raymond/Smith plasma 
(Raymond & Smith I1977t . the mekal model describing a 
Mewe/Kaastra plasma (e.g., Kaastra n992> and a simple power- 
law. The results of the spectral fits are presented in Tabled 

Because the jet is located deep inside the stellar body of 
the galaxy NGC 1052, we have to determine the amount of 
photoelectric absorption distributed along the line of sight. 
Using the thermal source models, we find an attenuating col- 
umn density value about one order of magnitude higher than 
the galactic foreground column density belonging to the Milky 
Way (Hartmann & Burton 1997 1. We attribute this additional 
X-ray attenuation to weakly ionized gas located inside the 
galaxy itself. The unabsorbed flux of the Raymond/Smith 
(Mewe/Kaastra) plasma in the range between 0.3 and 8keV 
corresponds to an intrinsic luminosity of L = 2.4(2.1) x 
lO"*" erg s"' at the distance of the source of 22.6 Mpc. 

Both plasma models fall below the measured count rate 
above 2.5 keV (see Fig.O suggesting some contribution of an 
additional hard X-ray component to the spectrum of the ex- 
tended emission. Due to the low photon statistics and the large 
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Fig. 4. HST structure map of NGC 1052. 
Dark regions represent dust obscuration 
while bright regions are locations of en- 
hanced emission. Overlaid is the pure nat- 
urally weighted MERLIN map showing the 
core of NGC 1052 with sub-arcsecond res- 
olution. Contours at (-1, 1, 2, 4, 8, 16, 32, 
64, 128, 256, 512, 1024) x 0.5mJy/beam 
are shown. The inlaid panel shows an en- 
larged view of the marked box around the 
area of the optical emission cone. 



Table 1. Best-fitting parameters for the spectral fits to the extended jet-associated emission 



Model<"> 


[10^2 cm-2] 


kT P"^ 
[keV] 


^abs 

[erg s"' cm"^] 


unabs 

[erg s"' cm""^] 


Xl, (^Vd.o.f.) 


(R&S)"''' 
(M&K)"''^ 


n 4Q+0.08 
"•""-0.02 


0.48!°:°^ - 
0.411°? - 

Q Q+0.6 


0.8x10^" 
0.7x10-'^ 
1.5x10-'^ 


3.9x10^'^ 
3.4x10^" 
2.0xlO-'3 


0.49 (10.8/22) 
0.53 (11.6/22) 
0.52 (11.5/22) 



" R&S: Raymond/Smith plasma. M&K: Mewe/Kaastra plasma, PL: Power-law; ^ Absorbing column density; 2 X 10^° cm ^ is 
set as a lower limit; ^ Photon index; ^ (Un)Absorbed X~ray flux between 0.2 keV and 8.0 keV 



uncertainties we have not tried to account for this additional 
spectral component by considering more complicated models. 

Because of the low number of counts from the jet, we can- 
not unequivocally dismiss the power-law model, which gave 
an acceptable fit to the X-ray data. However, we consider the 
power-law fit result of the X-ray jet as unreliable because of 
the unusually steep photon index"* (see Tablen*- 

4.2. The X-ray spectrum of the nucleus 

Below E - 2 keV a soft excess is detected in the nuclear spec- 
trum. Assuming that a fraction of the diffuse jet-associated 
emission originates within the immediate neighborhood of the 



The photon index T is defined as the power-law index of the spec- 
trum given in units of [photons s"' keV"']. It is related to the energy 
index a, used in radio astronomy, via: F = 1 - a 



nucleus, we used a hybrid model for the X-ray spectral approx- 
imation of the observed intensity distribution. 

/obs = (/diffuse + /nucleus X e--x'VH„o„,.) ^ e--xA'H„a,., 

where /diffuse is the X-ray spectrum of the diffuse jet-associated 
emission and /nucleus is the spectrum of the nucleus. Here, the 
photoelectric absorption produced by the torus is represented 
by e^^'^^mmu,, ^ while the X-ray absorption produced by the in- 
terstellar medium of the elliptical galaxy is e"°"><'*'ms»i > . 

We assumed a power-law type X-ray spectrum of the cen- 
tral X-ray source and modeled the soft component as a thermal 
plasma spectrum of the Raymond/Smith and Mewe/Kaastra 
type, (model I and II) respectively. We did not try to fit for 
the metallicities and fixed these during the spectral fitting pro- 
cedure given the poor statistics from the short integration time. 
Model la(IIa) assumes solar abundances of the Raymond/Smith 
(Mewe/Kaastra) plasma and allows the other parameters to be 
varied. In model Ib(IIb) the temperatures of the plasma and 
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Fig. 5. Top panel: X-ray spectrum of the extended X-ray jet emis- 
sion and the Raymond/Smith plasma model folded with the detec- 
tor response matrix (solid line). Middle panel: X-ray spectrum of the 
nuclear emission and the best fit of model lb (solid line, compare 
Table I4.2t . Bottom panel: The unfolded model lb spectrum repre- 
senting an absorbed Raymond/Smith plasma whose model parameters 
have been held fixed at the best fitting values to the extended X-ray 
emission spectrum (see Table0 plus a power-law with an additional 
absorber. 



the absorbing column density of the interstellar medium of the 
galaxy were fixed to the best fitting values to the extended X- 
ray emission (see Table 0. In model Ic(IIc) the metal abun- 
dances were additionally changed to 25% of the solar composi- 
tion. Such low metal abundances typically occur only in dwarf 
galaxies but are found in NGC 1052 from optical line measure- 



ments of Sil'chenko {1995} who reports values [Fe/H] 0.6 

for the nucleus and [Fe/H] 1 for the bulge of NGC 1052, 

corresponding to metallicities of about 25% and 10%, respec- 
tively. One way to accumulate a substantial amount of metal- 
poor material is a merger or close encounter with an extremely 
metal-poor dwarf galaxy. Evidence for such an event indeed 
exists from multiple observational approaches (see Forbes et 
al. 1200 H and references therein) in the case of NGC 1052. 

The statistically best fitting values are given in Table|2l For 
all six models we derive a luminosity of L = 1.4 x lO'*' erg s'^. 
Correcting for the effect of the two absorber model components 
increases the intrinsic (unabsorbed) value to L = (1-7 - 2.0) X 
10"*' erg s"' . The values of A^H(gai.) and kT for the plasma com- 
ponent in the models la and Ila are very similar to the values 
determined for the extended plasma emission. For these models 
we derive best fitting values for the absorbing column density 
of A^HCtorus) = (0.6 - 0.7) X 10^2 cm^^ as well as F = 0.2 - 0.3 
for the photon index. Similar values for both parameters result 
from fixing A^H(gai.) and kT as described above in the models lb 
and lib. Reducing the metal abundances of the plasma compo- 
nent in model Ic(IIc) forces the power-law to contribute more 
strongly to the soft part of the spectrum resulting in higher val- 
ues of A^H(toius) ~ 0.8 X 10^2 cm"2 and F ~ 0.3. These values 
are still rather low. Data from other X-ray observatories imply 
much higher values for A^H(toius) as well as for F (Weaver et al. 
1999 Guainazzi & AntonelH 1999 Guainazzi et al. l2000t . We 
note that the apparent discrepancy between previous X-ray ob- 
servations and the CHANDRA data concerning the photon in- 
dex might be due to the pile-up degradation (~ 30 %) of the lat- 
ter. The nuclear X-ray spectrum thus might appear artificially 
flattened. The determined absorbing column density, however, 
is not so sensitive to the pile-up effect, which mainly affects 
the hard part of the spectrum while the amount of absorption is 
determined in the soft X-ray regime. The statistical quality of 
all six fits is very similar, with values ofx^^^ ranging from 0.76 
to 0.84, and the low photon statistics (reflected in the values of 
^red < 1) do not allow us to dismiss any of these models. 

We note that a "patchy" absorber model in which a frac- 
tion of the central source is seen directly and only a part of 
the source is covered by the absorber might be a more realis- 
tic model for the nuclear X-ray emission of NGC 1052, given 
the results of Vermeulen et al. ('2003 1, who found that mov- 
ing VLBI components in both jets show complex light curves 
probably caused by substantial patchiness of the absorbing 
screen. Additionally, they find a complex Hi absorption line 
spectrum with compact clouds of absorbing material at differ- 
ent relative velocities and different locations along the parsec- 
scale jet structure. The absorbing column density derived from 
"patchy" absorber models is expected to be considerably higher 
(~ 2 X 10^2 cm"2) than the values derived from uniform ab- 
sorber models. However, given the low photon statistics, we 
omit the detailed discussion of such models with a larger num- 
ber of free parameters. 

5. Discussion 

The CHANDRA data provide for the first time direct evidence 
for jet-associated X-ray emission in NGC 1052. The diffuse, 
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Model*"' 


"H(gal.) 

[10^2 cm-2] 


kT 
[keV] 


2(0 


H(loms) 

[10^2 cm-2] 


p{e) 


abs 

[ergs"' cm"-] 


unabs 

[erg s"' cm"^] 


Xi, CrVd.o.f.) 


la 

lb 

Ic 

na 
nb 

lie 


0.49* 
0.49* 
51+0.09 

0.41* 
0.41* 


n yi /I +0 07 

^■^^ am 
0.48* 

0.48* 

0.28:°:«1 
0.41* 
0.41* 


1* 

1* 
0.25* 

1* 

1* 
0.25* 


r\ ^ A +0 S6 

0-74!o:39 
0-54!°:^« 
0.84!°:« 

0.82t« 


n m+0 10 

20+""' 

n. 90+0.09 

9^+0.06 
" -^-0.09 

27+°-'" 
"•^'-0.11 


2.3x10"'- 
2.3x10"'- 
2.3x10"'^ 

2.5x10"'^ 
2.3xlO"'2 
2.2x10"'- 


3.3x10"'^ 
2.9x10"'^ 
3.0x10"'^ 
3.4xl0"'2 
2.8xlO"'2 
2.8x10"'^ 


0.84 (15.9/19) 
0.76 (16.0/21) 
0.74 (15.8/21) 

0.89 (16.8/19) 
0.77 (16.1/21) 
0.78 (16.4/21) 



• Fixed value; » I: (R&S + PL*')'*^ II: (M&K + PL*0«'>», (R&S''"': Absorbed Raymond/Smith plasma, M&K'*<: Absorbed Mewe/Kaastra plasma, PL<*'>: 
(Absorbed) Power-law; ^ Column density of the diffuse absorber; ^ Metallicity (referring to the solar value); Column density of the compact absorber; ^ 
Photon index; ' Absorbed X-ray flux between 0.2 keV and 8.0 keV; 8 Unabsorbed X-ray flux between 0.2 keV and S.OkeV, i.e., Affj^^'^^, j and wj^Jl^^.^^j are set to 



extended X-ray emission can be best approximated with a ther- 
mal plasma model with kT ~ 0.4-0.5 keV. This temperature is 
consistent with the thermal component found earlier by Weaver 
et al. (1999) using ASCA and ROSAT daia. Its absorbed flux is 
only ~ 3% of flie nuclear X-ray emission but the intrinsic (ab- 
sorption corrected) extended emission might contribute up to 
14% to the total unabsorbed X-ray flux of NGC 1052. Because 
of the considerable pile-up degradation of the CHANDRA data, 
no firm conclusions on the photon index of the nucleus spec- 
trum can be deduced. The derived column density of hydro- 
gen towards the compact X-ray core (depending on the ap- 
plied model) of 0.5 - 0.8 x 10^^ cm"^ is in good agreement 
with the absorbing column density of ionized material towards 
the VLBI-jet derived by Kadler et al. (I2002> and Kadler et al. 
(in prep.). This suggests that the nuclear X-ray emission of 
NGC 1052 might be produced on the same scales as the parsec- 
scale structures imaged by VLBI at high frequencies. 

The detection of a diffuse region of X-ray emitting gas with 
a thermal spectrum and the same extent as the kiloparsec-scale 
radio jet suggests that jet-triggered shocks might play an impor- 
tant role in NGC 1052. In such a model the kinetic power of the 
radio jet is partially converted into X-ray emission. The optical 
morphology in the Hor filter substantiates this picture as was 
noted earlier by Allen et al. (1999 1. The alignment of the radio 
jet and the optical emission cone visible in Fig.0]implies that 
the ionization cone might be drilled out by the radio jets, result- 
ing in a predominantly shock-excited, conical narrow-line re- 
gion (see e.g., Dopita .2002j . Shocks might occur also on larger 
scales giving rise to the soft thermal X-ray emission associated 
with the radio jet/lobe structure in NGC 1052. A rough esti- 
mate (see Kadler et al. .2003 a) shows that the soft thermal X- 
ray spectrum associated with the radio jet of NGC 1052 can be 
explained in terms of the kinetic jet power being partially con- 
verted into X-ray emission originating in shocks driven into the 
ambient medium. (A more detailed model of the relation be- 
tween jet-driven shock-activity and the spectral shape of the 
extended X-ray emission in NGC 1052 will be discussed in 
a forthcoming paper) The comparison of the large-scale dis- 
tribution of radio emission in NGC 1052 between two epochs 
separated by ~ 15 years indeed shows activity on kiloparsec- 
scales. This substantiates the idea that shocks in the interstel- 
lar medium form at the working surfaces of active regions 
(hotspots and knots). Moreover, recent numerical simulations 



(e.g., Zanni et al. I2003> show that jets in radio galaxies can 
inflate over-pressured cocoons that drive shocks into the ambi- 
ent gas resulting in morphologies (in the case of weak shocks) 
very similar to what is observed in NGC 1052: a cavity of hot 
X-ray emitting gas in conjunction with a local deficit of X- 
ray emission around the hotspots. A deeper CHANDRA ob- 
servation with an improved photon statistic compared to the 
observation discussed here would provide both a higher sen- 
sitivity to the weak diffuse emission and a higher resolution. 
Additionally, the full resolution of CHANDRA of ~ 0.5 arcsec 
would allow one to study in more detail the connection between 
the knots in the diffuse X-ray emission and the optical emission 
knots. 

Acknowledgements. We thank G. V. Bieknell for helpful discussions 
and important suggestions. This research used the CHANDRA Data 
Archive (CDA) which is part of the CHANDRA X-ray Observatory 
Science Center (CXC) which is operated for NASA by the 
Smithsonian Astrophysical Observatory. We used CHANDRA data 
from an experiment planned and scheduled by G. P. Garmire. We 
made use of the data archive at the Space Telescope Science Institute 
which is operated by the Association of Universities for Research 
in Astronomy, Inc. under NASA contract NAS 5-26555. MERLIN 
is a National Facility operated by the University of Manchester at 
Jodrell Bank Observatory on behalf of PPARC. We made use of 
public MERLIN data from an experiment planned and scheduled by 
A. Pedlar. This research has made use of NASA's Astrophysics Data 
System. 

References 

Allen, M. G., Koratkar, A. R, Dopita, M. A. 1999, BAAS, 31, 893 
Dopita, M. A. 2002, Rev. Mex. Astron. Astrofis., 13, 177 
Falcke, H., Wilson, A. S., Simpson, C. 1998, ApJ, 502, 199 
Forbes, D. A., Georgakakis, A. E., Brodie, J. R 2001, MNRAS, 325, 
1431 

Gabel, J. R., Bruhweiler, F. C., Crenshaw, D. M., Kraemer, S. B., 

Miskey, C. L. 2000, ApJ, 532, 883 
Guainazzi, M., Antonelli, L. A. 1999, MNRAS, 304, L15 
Guainazzi, M., Oosterbroek, T., Antonelli, L. A., Matt, G. 2000, A&A, 

364, L80 

Hartmann, D., Burton, W. B. 1997, 'An Atlas of galactic neutral hy- 
drogen", Cambridge (New York: Cambridge University Press) 
Heckman, T. M. 1980, A&A, 87, 152 

Kadler, M., Ros, E., Kerp, J., et al. 2002, In Proceedings of the 6th 
European VLBI Network Symposium, Ros, E., Porcas, R. W., 



8 



M. Kadler et al: Jet emission in NGC 1052 at radio, optical, and X-ray frequencies 



Lobanov, A. P., Zensus, J. A. (eds.), (Bonn, Germany: MPlfR), 
p. 167 

Kadler, M., Ros, E., Kerp, et al. 2003a, New Astronomy Reviews, 47, 
569 

Kadler, M., Ros, E., Zensus, J. A., Lobanov, A. R, Falcke, H. 2003b, 
SRT: the impact of large antennas on radio astronomy and space 
science, N. D'Amico, F. Fusi Pecci, I. Porceddu, G. Tofani (eds.), 
(Bologna: Societa ItaUana di Fisica), p. 219 

Kameno, S., Sawada-Satoh, S., Inoue, M., Shen, Z., Wajima, K. 2001, 
PAST, 53, 169 

Kaastra, J. S. 1992, An X-ray Spectral Code for Optically Thin 

Plasmas (Internal SRON-Leiden Report, updated version 2.0) 
Kellermann, K. 1., Vermeulen, R. C., Cohen, M. H., Zensus, J. A. 

1999, BAAS, 31,856 
Knapp, G. R., Faber, S. M., Gallagher, J. S. 1978, AJ, 83, 139 
Krist, J., Hook, R. 1999, The Tiny Tim User's Guide, Ver. 5.0 

(Baltimore: STScI) 
Malizia A., Bassani L., Stephen J. B., Malaguti G., Palumbo G. C. C. 

1997, ApJS, 113, 311 
Mc Dowell, J. C. 1994, in Einstein Obs. Unscreened IPC Data Archive 
Pogge, R. W., Maoz, D., Ho, L. C, Eracleous, M. 2000, AJ, 532, 323 
Pogge, R. W., Martini, R 2002, ApJ, 569, 624 
Raymond, J. C, Smith, B. W. 1977, ApJS, 35, 419 
Risaliti G., Elvis M., Nicastro F. 2002, ApJ, 571, 234 
Sil'chenko, O. K. 1995, Astron. Letters, 21, 283 
Sugai, H., Malkan, M. A. 2000, ApJ, 529, 219 
van Dokkum, R G. 2001, PASR 113, 1420 

Vermeulen, R. C, Ros, E., Kellermann, K. I., et al. 2003, A&A, 401, 
113 

Weaver, K. A., Wilson, A. S., Henkel, C, Braatz, J. A. 1999, ApJ, 520, 
130 

Wrobel, J. M. 1984, ApJ, 284, 531 

Zanni, C, Bodo, G., Rossi, R, et al. 2003, A&A, 402, 949 



